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Abstract 

We evaluate the corrections to the mean field values of the x and 
z exponents at the first order in the e-expansion, for T = T c . We find 
that both x and z are decreasing when the space dimension decreases. 

We want to investigate the purely relaxational dynamics of a short-range 
spin glass model for T — > T 6 + , in the framework of the e-expansion. The 
dynamical properties of the model in the mean field theory are very different 
from those of the models whose dynamical properties are usually investigated 
in the literature. These new feature make the computation of the dynamical 
critical exponents much more involved than that of the usual. 

In a previous work, ||, we evaluated the Gaussian dynamical fluctuations 
of the order parameter around the MF limit. The aim of this letter is to 
pursue this analysis, by considering the 1-loop correction to the mean field 
(MF) theory in a renormalization group calculation. At the first order in e, 
unlike Zippelius ||, we find results that disagree with the conventional Van 
Hove theory, because we obtain a correction to the kinetic coefficient already 
to the lowest order in the loop expansion. In this work, first we evaluate the 
correction to MF value of the critical exponents x, that describes the critical 
slowing down of the dynamical order parameter at the critical point, then we 
evaluate the correction to the z exponent that describes the critical slowing 
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down of the dynamical spin glass susceptibility Xsg- Finally, we check that 
the scaling low, which connect these two exponents, is verified. 

We study the soft-spin version of the EA model given by the Hamiltonian: 



(ij) i i 



(i) 



where are random Gaussian interactions between the nearest-neighbors 
sites, with zero average and mean square fluctuations [(Jij) 2 ] = j 2 /n (n is 
the coordination number). Purely relaxation dynamics is introduced by the 
Langevin equation: 



.^(t) d(/3H) 



dt 



ds t {t) 
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where r\ is the usual Gaussian noise with zero average and variance (r)i(t)r)j(t')) = 
Y^5ij5(t — t'). The interesting physical quantities in MF theory are the aver- 
aged response and correlation functions, defined respectively: 

d(si(t)} v 



G{t - 1') 



dhi{t>) 



C(t-t') = [{s i (t)s i (t')) r ,} 



(3) 
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where the angular brackets (..) v refer to averages over the noise and the 
square brackets over quenched disorder. 

Moreover, considering the Gaussian fluctuation, we can define the dy- 
namical spin glass susceptibility as follow: 



Xsc{i - j;h - t u t 2 - U) 



d(s i (t 3 )) ri d{s i (t 2 )} 1 
dhj(ti) dhj(t A ) 
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The dynamical scaling implies that the decay of C(t) is governed by a 
characteristic time r, which diverges at T c , as, for long t, we can write: 



C(t) 



(6) 



where q + is the universal scaling function in the region T — > T c + . The relax- 
ation time divergence, at the critical point, is connected to the correlation 
length divergence through the dynamical exponent z: 



(7) 
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According to the scaling hypothesis, this exponent is related to the slowing 
down of the spin glass susceptibility, and we expect to have: 

Xs G (k,u)=u*?f(k*/u) (8) 

The MF behaviour of this model (n = N, number of spins, long-range 
limit), in the critical region, is well known 0, ||. In the low- frequency 
limit the response and correlation functions are respectively: 

G(u) = (l - V=i^) , (9) 
C(u) = 2 (10) 
while, in the Gaussian approximation, the spin glass susceptibility is 

Xsc{k, wi, uj 2 ) = 2 /-^-, /—■ — ( n ) 

The MF value of the dynamical critical exponents x and z, as known, are 
1/2 and 4 respectively 

To deal with Langevin disordered dynamic theory as usual we use the 
dynamic functional integral method 0. In this formalism, it is conventional 
to introduce an auxiliary field Sj(t) and to define an effective Lagrangian of 
an Hubbard-Stratanovich field Q° (t,t'), @, ||, such as: 

2(Q a l o (t 1 ,t 2 )) L{Qa0) = M(ti)tf(t2)) v }j (12) 

where the two component vector field is defined: 

(f>?=(isi,8i). (13) 

In considering the correction to MF approximation we derive the following 
Lagrangian, as series expansion around the order parameter saddle-point 
value Q^(M'), I, EH, & 

L{8Q^) = - Y,Y,K- J 1 5Qf(t 1 ,t 2 )A^ 5 5Q?(t 1 ,t 2 )+ 

\ E Y:SQf(t 1 ,t 2 )C a ^\t 1 ,t 2 ,t 3 ,U)SQf(t 3 ,U)+ 

h E E C a ^(h, t 2 , t 3 , u, t s , t e )8Q?(t u t 2 )5qfr s {h, U)5Q?{t s , t 6 ) . 

(14) 
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For the structure and the meaning of each term of (|I4"D, the reader is referred 
to 0. The non local connected propagators of the theory are defined as 
follow: 



[<^(«l)rf(* 2 )>,]j[W(*3)^(«4)>,]j 

= 4 Ei^h T,(K~ 1 )*(8Q?(ti, h)5Ql\t 3 , U)) L{QaB) 

l k 

-2(K- 1 ) ij A afhfS 5(l - 3)5(2 - 4) . (15) 



In we evaluated the critical behaviour of these propagators for each 
combination of the indices a,/3, 7, S and in any time interval, when the 
cubic interactions is vanishing. We write down the general structure of 
G 2211 (u)i, u 2 , <^3, ^4) which is present in several of the following 1-loop func- 
tions: 

G 2211 (fc;Co>i,Co>2,^3,^4) = [S(ui + LU 3 )5(uJ 2 +W4) + 5(lU 2 +Co> 3 )(5(cl>i + OJ4)] 

x2211 

G 2211 (k; uiu 2 ) + 5(ui + to 2 + to 3 + uj 4 )G (k;ux,U2iU$,U4) (16) 

where: 

G 2211 ^,^) = 1 , (17) 

x2211 . _ _, > 

G (k; ut, u 2 , u 3 , u A ) = [C(ux)G 2211 (k] -wi, w 2 ) + CV 2 )G 2211 (/c; £*>i, -ofe)) x 

G 2211 (k; tvx, u 2 )G 2211 (k- -u 3 , -uj 4 )g r (k; u x + u 2 ) (18) 
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g r (k] ujx + lo 2 ) = (ck 2 + yf-ifa + u^ Fx f ^- ^^172 J i ( 19 ) 

where -Fi(x) is an homogenous function of the We recall that in this 
formalism the time-dependent spin glass susceptibility (||) is: 

Xsc(i ~ j]t 3 -tx,t 2 -t 4 ) = [(s i (t 3 )r] j (ti)) v (s j (t 2 )r] i (U)) ri ] j = G 1221 (i - j ; t 3 -tx,t 2 -t 4 ). 

(20) 

Let us consider the 1-loop correction to the "free" theory. We intend to 
use the propagators derived in || to evaluate the contribution of the 1-loop 
Feynman diagrams to the mean value Q a @, to the bare propagators G a ^ 6 ) 
and to the bare cubic vertices. 
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Figure 1: 



We consider the 1-loop Feynman diagrams as a g-series expansion, by 
using the correspondent propagators. In analogy with the static, we can 
guess that the (^-dependent part of the 1-loop corrections is not singular at 
the critical point, as soon as g ^ 0. The physical quantities have to be not 
affected by the value of g, provided that it is not zero. The behaviour at 
g = is rather different and the e-expansion starts from D — 8, fl2 |. 



Concerning the 1-point function, let us consider the response function, 
that as a consequence of eq. (0) is: 

G(u) = 2<g£») (21) 

The diagrams that occur in 1-loop correction to the MF value of the x ex- 
ponent are shown in Fig. [1]: two continuous lines represent a bare factorized 
in time propagator 

G aM& {k- Wi, W 2 ) [6(U! + ^3)5(^2 + w 4 ) + 5{u x + u 4 )5(lu 2 + CU 3 )} , 
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Figure 2: 



three lines the bare connected in time propagator 

G (k;wi,w 2 ,u s ,u 4 )5(uj 1 + u 2 + u 3 + u; 4 ) 

and, finally, the triangle in the center of the diagrams is a factorized in time 
cubic vertex. We succeed in evaluating the singular behaviour of the diagrams 
A\ and A 2 , while we have to use a trick to take into account the contribution 
from the others. By using the series expansion in g of the propagators in the 
loop, see H, we obtain a correspondent g series of the two 1P.I. diagrams, 
but the zero term is missing. Let us add and subtract the term we need. As 
in the static case (where we deal with just the momentum variable k, and we 
easily manage to resum the ^-series), we can suppose that the resumation of 
the series give a non singular behaviour at the critical point (the presence of 
g remove the pole at zero momentum for T = T c ) . We left with the diagrams 
Ai, A 2 and with the two zero-order in g diagrams that we need for the series 
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resumation (with negative sign), (see Fig. [2]): 



d d k 



(27r) d c k 2 + 2VkJ 



d d k 



(2vr) d ck 2 + Viu 



(22) 



The propagators involved in these "new" diagrams are represented, in Fig. [2], 
with two crossed lines. By evaluating the previous integrals for d = 6, we 
find that, at the first order in the loop expansion, the response function is 



G(u>) 



l - (-iuY 
l 



G(lo) + 2u (5Q 2 \lu)) 



V- 



■IUO 



(2vr) ! 



I(uj = 0) + -y/-iu \n(u) 



(23) 



As for the statics u has been introduced as expansion parameter in the cubic 
vertices, to apply the renormalization group method for critical phenomena. 
The factor u 2 in the equation (|23|) is due to the fact that the saddle point 
response function @ is of order (-). 

In the same way, we derive the flux equation of the IP. I. vertex functions 
which allow us to determine the fixed point below 6 dimension. As for the 
conventional Langevin dynamical theories, the IR stable fixed point below 
6 dimension is the same that the static one (i.e., we find the same relevant 
diagrams) : 

, _ (2tt) 6 £ 
k. 2 



(24) 



The series expansion in e for the static exponents is evaluated up the 3th 
order in Ref. 
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Let us determine the first order correction for the dynamical ones. Sub- 
stituting the fixed point value (p4|) into fl2"3|) we obtain for x: 



(25) 



1 / e 
x = - 1- - 

2 V 4 



In the same way, we find that the relevant 1-loop IP. I. diagrams for the Self- 
Energy S(a; 1 ,u; 2 ), shown in Fig [3], give the following contribution to the 
spin glass susceptibility 



1 



1 



u 2 f d d k ( 

2 J (2ll) d \ C(p - k) 2 + 2y/lUhck 2 + y/iul + y/^EJ^ 



(26) 
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B 



u 



d d k 



1 



1 



C : -2 — 



2 J (2ir) d \c(p - k) 2 + 2^vJ 2 ck 2 + y/iu^ + ^f^vJ 2 , 
d d k ( 1 1 \ 



2 7 (2vr) d \c(p - k) 2 + yfiul ck 2 + 



D : -2 



1 



£:^2 



(2ir) d \ c(p — k) 2 + \/—iuJ2 ck 2 + 
d d k 1 



2 q» 2 + (2n) d ck 2 + 2^/Iuj^ 

u 2 1 /■ d d fc 1 

2 cp 2 + 2y/iuo[J (2ix) d ck 2 + ^fw[ 

,2 



G 



1 



d d k 



1 



2 qo 2 + 2y^i^J {2it) d ck 2 + 2yf=i^~ 2 



u 2 
-2— 2 



d d k 



2 cp 2 + 2sf^Zr 2 J {2Ti) d ck 2 + t/=iu~ 2 



(27) 
(28) 
(29) 
(30) 
(31) 
(32) 
(33) 



At the first order in e, the sum of the contributions to z exponent of the 
diagrams A, B, C, D is vanishing. On the other hand, from the sum of the 
diagrams E, F, G and H we obtain the following contribution to the 1P.I. 
2-point function: 



p=0 

u 2 k 



4 (2tt) 



(y^iuJi \n{u x ) + \f-iuj 2 In (^2)) (34) 



At the fixed point u* given from (0), we obtain the following correction 
to z exponent: 



(2-77) 
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(35) 



The scaling relation between the exponents x and z, that we recall to be 

d-2 + r] 



x 



2z 



(36) 



is verified, at the first order in e. 

Numerical simulations for the exponents x and z can be found in the 
literature, the values are z ~ 7 and x ~ .06 in dimensions -D = 3, [13], and 
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z rs 5 and x ~ .15 in dimensions D = 4, [I3|. Our prediction states that 
both the values of x and z are decreasing when the dimension decreases. 
This is true for x but not for z. The apparent discrepancy that we have with 
the behaviour of z should not worry us. In fact, also in the static case, the 
critical exponents for spin glasses have a badly convergent e-expansion and 
the prediction of this expansion can be hardly applied in 3 or 4 dimensions. 

A numerical study of what happens in 5 dimensions is necessary. More- 
over it should be noticed that usual arguments imply that our computation 
predicts, without ambiguities, that the logarithmic corrections in 6 dimen- 
sions are such to decrease the effective value of the exponents. 
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CAPTIONS FOR ILLUSTRATIONS 



Fig.l 1-loop diagrams that contribute to the x exponent. 

Fig. 2 g-series resumation of the 1-loop diagram contribution to the x exponent. 
Fig. 3 g-series resumation of the 1-loop diagram contribution to the Self-Energy. 
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